
NUMERICAL ANALYSIS OF 
STIFFENED SHELLS OF REVOLUTION 

Volume VI of VII 


by V. Svalbonas and P. Ogilvie 


Prepared by 

GRUMMAN AEROSPACE CORPORATION 

Bethpage, N.Y. 11714 

for George C. Marshall Space Flight Center 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • SEPTEMBER 1973 



TECHNICAL REPORT STANDARD TITLE PAGE 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA CR-2273 


4. title ano subtitle 

Numerical Analysis of Stiffened Shells of Revolution 
Volume VI of VII * 


7. AUTHOR! S) 

V. Svalbonas and P. Ogilvie 


9. PERFORMING ORGANIZATION NAME ANO ADDRESS 

Grumman Aerospace Corporation 
Bethpage, New York 11714 


12. SPONSORING AGENCY NAME AND AODRESS 


3. RECIPIENT'S CATALOG NO. 


5. REPORT DATE 

September 1975 


6. PERFORMING ORGANIZATION CODE 


6. PERFORMING ORGANIZATION REPORT » 

M381 


llO. WORK UNIT NO. 


11. CONTRACT OR GRANT NO. 

NAS 8-27468 

13. TYPE OF REPORT & PERIOD COVERED 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


Contractor Report 


14. SPONSORING AGENCY CODE 


IS. SUPPLEMENTARY NOTES 

•fa . 

A list of all seven volumes, by sub-title, is on the back side of this page 


16. ABSTRACT 

This manual contains engineering programming information for the STARS-2V 
(Shell Theory Automated for Rotational Structures -2V (Vibration)) digital computer 
program. The report is written for the engineer who will need to make small 
alterations to the program, such as incorporating a new geometry, or altering a 
table size, to fit his specific needs. The sections of the manual each cover one 
major subroutine. These sections are further subdivided in the following manner 
where applicable: 

A. Subroutine description. 

B. A list of pertinent engineering symbols and their FORTRAN coded 
counterparts. 

C. Subroutine flow chart. 

' h ; ■ 1 

D. Subroutine FORTRAN listing. 



19 . SECURITY CLASS1F. (of thU report 

20. SECURITY CLASS1F. (of thU page) 

21. NO. OF PAGES 

Unclassified 

Unclassified 

192 


\1SFC- Form 3 29 2 (Rev December 1911) 


For sale by National Technical Information Service, Springfield, Virginia 12151 





















VOLUME 


VOLUME 

VOLUME 

VOLUME 

VOLUME 

VOLUME 


I. Theory Manual for STARS-2S, 2B, 2V Digital Computer Programs 

II. Users' Manual for STARS-2S - Shell Theory Automated for 

Rotational Structures - 2 (Statics), Digital Computer Program 

III. Users! Manual for STARS-2B, 2V - Shell Theory Automated 

for Rotational Structures - 2 (Buckling, Vibrations), Digital ' 
Computer Programs 

IV. Engineer's Program Manual for STARS-2S - Shell Theory 
Automated for Rotational Structures - 2 (Statics), Digital 
Computer Program 

V. Engineer's Program Manual for STARS-2B - Shell Theory 
Automated for Rotational Structures -2 (Buckling), Digital 
Computer Program 

VI. Engineer's Program Manual for STARS-2.V - Shell Theory 
Automated for Rotational Structures -2 (Vibration), Digital 
Computer Program 


VOLUME VII. Satellite Programs for the STARS System 



ACKNOWLEDGEMENT 


The authors wish to acknowledge the overall aid of Mr. J. Key, especially 
in the debugging phases of the UNIVAC program version. 



INTRODUCTION 


This manual presents a general description of the STARS-2V digital computer 
program. FORTRAN IV is used exclusively in writing the various subroutines. 
The- execution of this program requires the use of eleven temporary storage 
units . 

i' 

The program was initially written and debugged on the IBM 370-165 computer 
and then converted to the UNIVAC 1108 computer, where it utilizes approxi- 
mately V7>000 words of core. Only basic FORTRAN Library routines are 
required by the program, these being: sine, cosine, absolute value, 
square root, min 0, max 0, amax 1, amin 1, and random. 

For ease and speed in usage, the Table of Contents on the following page 
has also been laid out to present the call sequence of the program. 
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SUBROUTINE MAIN 

MAIN is .'the control link for the entire program. Sizing values are read 
into the program, as well as information for eigenvalue or determinant 
accuracy' calculations, assumed frequencies and resonance factors, and'^he 
material property tables. Calls are made to subroutines RIEMAN and SEGMAT 
once for each segment in a region; then subroutine REGMAT is called. This 
procedure is executed once for every region in the structure. Finally a 
call to subroutine STRMAT is made. If the analysis involves prestressing, 
subroutines INITAL and LEBEGE are' called. Otherwise the control link 
loops back and resets the clues for the calculation of dynamic stiffness 
matrices. 

In the calculation loop for the dynamic stiffness matrices the program 
proceeds as above to STRMAT. Next EIGVAL or DETERM- are called to calculate 
eigenvalues or evaluate a determinant as per user input. If these 
calculations satisfy convergence criteria, subroutines BCVECT, INITAL and 
LEBEGE are called. If convergence criteria are not satisfied the program 
will loop until such satisfaction is obtained. 

In a multi-harmonic eigenvalue search the whole procedure loops on the 
harmonic number. 

several counters in this control link. These are defined, as 

Counts the calls to subroutines RIEMAN and SEGMAT, from 
1 to the number of segments within a region. 

Counts the calls to subroutine REGMAT, from 1 to the 
number of regions in the structure. 

Cycle counter (counts harmonics). 

Pass counter (counts passes necessary to converge in a 
single harmonic). 

and overlay listings are included in this section. 


There are also 
follows : 

NSC 

NRC 

ICYC - 
NPASS - 

The block data 



MAIN 
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SUBROUTINE RIEMAN 


This subroutine link assembles the data tables for use in the integration 
procedure. The program has the capability of handling two loading conditions 
(due to rotation or static prestress). Temperature load can be included 
in the static prestress. 

The subprogram link, RIEMAN, utilizes the subroutines SETUP, ROBOT, DIFl r 
or DIEF2, to integrate the differential equations of each segment independently, 
under arbitrary load conditions. The results of the integrations of each 
segment are stored in the YCORR array in RIEMAN, and represent the stiffness 
and deflection coefficients of each segment. 
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SUBROUTINE SETUP 


SETUP is a' double entry subroutine called from RIEMAN. It is a mixed precision, 
numerical integration routine, with automatic selection of a variable integration 
step size, which utilizes fifth order Runge-Kutta equations to obtain the solution 
for first order differential equations. 

SUBROUTINE MAGIC 


MAGIC is an alternate entry point to subroutine SETUP. 




CALCULATE K th 
RUNGE-KUTTA FORMULA 









SUBROUTINE ROBOT 


This subroutine is used by RIEMAN to calculate geometric and . load 
coefficients for use in the differential equations. With reference to 
geometry, all the necessary radii are calculated, as 'well as the stiffness 
coefficients of the various shell wall constructions. Thermal load moments 
and direct forces are also calculated" from direct temperature input. Loads 
due to rotation, and the shell wall areal mass are also calculated by ROBOT. 

All the above values are passed back via the label common area EQUAZN. 

In the case of a special point input geometry the ROBOT routine calls 
GEOMET . 

Subroutines GEOMET, PLIGO, PLINE 

Starting from a set of z, r points these subroutines calculate the necessary 
radii of the shell curves using spline fits. 
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FORTRAN CODE 

ENGINEERING SYMBOLS (REF. 1) 

RO 

r o 

R1 

r i 

R1DOT 

It a 


1 ><t> .. . 

. CS :• > 

. ; COS <p 

SN 

■ sin <j> 

. A 

a 

C ■ 

■ c. 

XN 

, . n 

F2 

f 2 

F3 ' ' ‘ 

’ f 3 ' •' '' ' “ 

TAN; TN 

tan <j) 

SEC 

sec <j> 

TII 

T. . 


n 

TIK 

T. .. . . 


• 1C . . •_ 

TOK 

T 


oc 

TOO 

T , 


oo - 

TEFREE 

T 

HI 

h. 


l 

HO 

h 


o 

T 

t 

TI 

t. 


l 

TO 

t 


o 

SNSQ 

• 

sm <p 

CSSQ 

2 

COS <p 

CN 

, cos <J> sin <p 

X1CS 

l/cos <t> 

X1SN 

l/sin (j> 

R2 

r 2 


BETA 
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FORTRAN CODE 


ENGINEERING SYMBOLS (REF. 1 ) 


X1ROSN 

X1ROCS 

CSX1RO 

CSX1R1 

CSXIR2 

SNX1RO 

SNXIR1 

X1R1 

X1R2 

X1R1SQ 

X1ROSQ 



1/ r Q sin 4> 
l/r Q cos <j> 
cos 4>/r Q 
cos cJj / r^ 
cos 4 >/r 2 
sin <(> / r 
sin ij>/r 




2 
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ROBOT 


ENTRY 
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SUBROUTINES DIF1 AND DIEF2 


These subroutines are called in RIEMAN as necessary. DIF1 contains the 
differential equations for the THIC and ST clues , while DIFF2 contains the 
differential equations for the RWA and ISG clues. Geometry clues, 
trigonometric values and predicted values of the differential equation 
variables are passed via label common area, EQUAZN, to subroutines DIF1 
or DTFF2. The coefficients for the mass and load terms, XI, X2, X3, and 
K are identified depending upon the pass number and convergence criterion. 

The specific derivative equations and auxiliary equations are contained in 
these subroutines. The values of each derivative equation, YDOT, and each 
auxiliary equation, YA , are returned to RIEMAN via label common EQUAZN . 

A special equation counter, I, is used in these subroutines, which counts 
in increments of eight. The first eight values of I, 1 through 57 (in 
increments of eight), correspond to the eight sets of initial conditions 
required to compute the segment stiffness matrices in subroutine SEGMAT. 

The subsequent values of I, 65 and 73 maximum (again in increments of eight) 
correspond to the computation of each set of eight equations for each 
loading condition (two conditions corresponding to rotation and static 
prestress respectively) . 
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FORTRAN CODE 


XN 


ENGINEERING SYMBOLS 


n 


(REF. 


1) 


YDOT (I) 

YDOT (I + 1) 

YDOT .(I + 2) 

YDOT (I + 3) 
YDOT (I + 4) 

: YDOT (I + 5) 
YDOT (1+6) 

YDOT (1+7) 

YPRED (I) 

YPRED (I + l) 

YPRED (I + 2) 

YPRED (1+3) 

YPRED (1+4) 
YPRED (1+5) 
YPRED (I + 6) 
YPRED (1+7) 

YAMPT 

YANTH 

YAMTH 


T 

<pd ,<p 

dT *e 

ds 

n <m' 

dN 

i 

ds 

J <j) ,<j) 

dJ A 

2 

ds 

M. A 

<P ><P 

- dM 

2 

ds 

U,. 

<P 

dU 

ds 


d V 

* » X 

ds 

% 

dW 

ds 

V* 

% 

ds 


<f>0 


N, 


a i 

U 

V 

w 


M 


<f>0 




/ 
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ENGINEERING SYMBOLS (REF. 1) 



Non-Linear Redefinitions (Ref . 1 ) 


YDOT (I+2-) 

i 

*J . . 

<t> ><f> 

) 

d*J 

i 

ds 

YPRED (1+2) 

*JT" 

<P 


YAJPH . 

J A 1 
<t> 1 


XNL 

a 


XNPHI 

\ 

1 


) 
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DIF 1 


ENTRY 
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SUBROUTINE SEGMAT 

The results of the subroutine link, RIEMAN, are passed through the label 
common area, LfCORR, to this subroutine. SEGMAT places the elements of 
the YCORR array into several double-subscripted arrays, forms some 
coordinate transformation arrays, and calls subroutine SREVN2 for double 
precision matrix inversion. 

As a result of appropriate matrix operations this subroutine produces a 
segment stiffness matrix, the XKS array, and a segment load matrix, the 
XIS array, for each segment. In passes other than the first pass, first 
cycle, the XLS array is not calculated. SEGMAT also orients -each segment 
into the global coordinate system of the structure as a result of the 
matrix operations. 

Subroutine SREVN2 

SREVN2 is a subroutine called by SEGMAT to invert a real, double -precis ion, 
in-core matrix utilizing Gauss-Jordan elimination with partial pivoting. 
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FORTRAN CODE 

ENGINEERING 

SYMBOLS (.REF. 1 ) 


SNI 

si 




SNJ 

sj 




CSI 

ci 




CSJ '* 

cj 




A' MATRIX 


' i ft 

I 0 
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j JFT 



Bi MATRIX 
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1 A 
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; x i 
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X 3. 
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j 0 

j 0 
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1 Ip . 
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j 0 
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XKT MATRIX 
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XMAX MATRIX 
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XKS MATRIX 
XLS MATRIX 
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SUBROUTINE REGMAT 


The segment stiffness matrices, XKS,' and the segment load matrices, XL3 , are 
passed from SEGMAT to REGMAT via Tapes #2 and and are placed in the XKRTOT 
array and the XLRTOT array, respectively. If kinematic links occur between 
segments in the region, the XKRTOT array and the XLRTOT array are modified 
to represent the situation. In the case of discrete rings the routine RINGER 
is ’called and provides the necessary matrices. 

A horizontal and vertical partitioning of the XKRTOT array occurs while the 
XLRTOT array is subjected to a horizontal partitioning only. Appropriate 
matrix operations are performed upon the partitions of each array, thus 
reducing the size of the region stiffness and load matrices and resulting 
in increased program capacity. The results of these manipulations are the 
region stiffness matrix, XKR, and the region load matrix, XLR. In passes 
other than the first pass, first cycle, the load calculations are not performed 

Subroutines Called from REGMAT 

Subroutine SWITCH ; Is a routine used to arrange a matrix in a form convenient 
for use by another routine employing a positive definite method for solving 
linear algebraic equations. 

Subroutine CHASE : Is a routine used to obtain the solution X of the linear 

system AX = Y, given at least one right side of Y and the positive, definite, 
symmetric, real coefficient matrix A. 

Subroutine FUTILE: Is a routine called from CHASE and used to obtain the 

factorization of the positive definite, real, symmetric matrix A into the 
product of a lower triangular matrix and its transpose by utilizing a 
Cholesky decomposition. 

Subroutine TRIEQ : Is a routine called by CHASE or EIGVAL to solve a triangular 

system of algebraic equations. 
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FORTRAN CODE 

ENGINEERING SYMBOLS (REF. 

SKL MATRIX 

1 

[sklJ 

SKLTR MATRIX 


[skl] t 

XKRTOT MATRIX 





K, 21 j K, 22 


XLRTOT MATRIX 


SKL22 MATRIX 


REGTOT MATRIX 


STORE MATRIX 


XK11 PARTITION 


XK12 PARTITION 


XK22 PARTITION 


XK21 PARTITION 
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FORTRAN CODE 


ENGINEERING SYMBOLS (REF. . 1 ) 


XL1 PARTITION 


XL2 PARTITION 


XKR MATRIX 


N 


XLR MATRIX 
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SUBROUTINE RINGER 


This subroutine reads the discrete ring geometric data, and temperatures, 
and forms the ring stiffness, mass, and thermal load matrices. These 
matrices are passed back to either of subroutines REGMAT or.STRMAT (see 
next) as necessary, for incorporation into the region or structure matrices, 
respectively. The effects of static prestress are considered, and the load 
calculations are omitted in passes involving eigenvalue extractions. 

The calculations in RINGER account for the eccentricity of the ring centroid 
from the base shell wall, and the offset of the ring centroid from the 
shear center. 
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FORTRAN CODE 

ENGINEERING SYtCBLS (REF. l) 
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SUBROUTINE STRMAT 


The region stiffness matrices, XKR, and the region load matrices, XLR, are 
passed from KEGMAT to STRMAT via Tape #4 and Tape # 8 , and are placed in 
the XKSTOT array and the XLSTOT array, respectively. A matrix, BCD, is 
formed to represent the boundary conditions, and, if kinematic links 
occur between regions, the RKL matrix is developed to represent this 
situation. The subroutine RINGER is again called for discrete ring 
matrices. In passes other than the first pass, first cycle, the load 
calculations are not performed. 

As a result of appropriate matrix operations, a reduced structure stiffness 
matrix is formed. .In passes other than the first pass, first cycle, STRMAT 
terminates here, in the first pass, first cycle, the following calculations 
are made. Subroutine FLEX, a routine identical to SREVN2 (except for being 
in single precision) with the name changed due to the structure of the 
OVERLAY option, is called to invert the structure stiffness matrix thus 
producing A, the flexibility matrix for the structure. The region end 
deflection array, DRE, is produced as the result of another set of matrix 
operations. 
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FORTRAN CODE 

ENGINEERING SYMBOLS (REF. l) 

BCD MATRIX 

[ ) 

BCT MATRIX 

[BC] T 


A 

XST MATRIX 

( ]. ' 



XKF MATRIX 

Wf 


A 

A MATRIX 

[a] f , 


A 

XSL MATRIX 

f l] 


..... .. . - L J ip 


A 

. XLS ARRAY 

- (If - 


A 

XFL ARRAY 

1 f 1 f 

DRE ARRAY 

- f A } 


l it 
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STRMAT 
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RETURN 1 504120 

END 1504130 





END 1600491 



SUBROUTINE INITAL 


As a result of the matrix operations performed in REGMAT, the SKL22, the 
XK2221, and the XK22I 2 arrays for each region are passed to INITAL. The 
XK1112 and XL1 arrays for each segment, resulting from the matrix procedures 
in SEgMAT, are also passed to INITAL. The region end deflection matrices, 
DRE, which were formed in STRMAT or BCVECT are transmitted to INITAL. 

Following appropriate matrix operations upon these arrays, the force initial 
conditions, the FICS array, and the deflections initial conditions, the 
DICS array, are produced. These arrays combine to form the YICS matrix, 
which contains the true initial conditions for the structure to be analyzed. 

The pertinent counters in the subroutine are: 

NS = segment counter 
NR = region counter 
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FORTRAN CODE 


XK2221 MATRIX 


ENGINEERING S YMBOLS 

V 

C 22 



(REF, l) 

A 


A ,-i 


XK22L2 MATRIX . 

1 ] 

DSE ARRAY 

w 

XK1112 MATRIX 

[“.* 


r i T 

ROTD MATRIX 

1 idtJ 

DICS ARRAY 

{6(D) 

XL1 ARRAY 

{i(D) 


r i T 

ROT MATRIX 

1 iftI 

FICS ARRAY 

{f(D) 
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DHTEMP = 0 


NR 
= t 


J 





FICS = 

ROT * TDSE 


i)SE = 

K 1 1 1 2 • DSE 


RESET TDSE - 
TDSE + XLt 




FORM ROT 


108 
























\ 


o o OOOOOOOOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO©*-* 
•Hf\ir^>4-u^or~cocro^cNjm-^iA>or^coo'0^rgro»rin>Oh-ooo'0»- < r\jf^'ti/''0^*a3cro*-«iNjf T ^»rir>'or^co(7'0»~ < (Mf<'-<7-i/>sor^cocra' 
0O00oo000HH^HH-irtHM-<(\(N(rg(MiNi(N(Mc\MV(N( i nmpn(fwnrnnnmm^^^-r' 1 r^^'j'^^uM/M/>inuMmPirmirit> 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOUOOCJOOOOOOOOOOOOOOOOOO 


to 

3 


x oc 

o • 


Q. tO 
50 Z> 

>- a 
to < 
UJ to 
QC * • 


— to 

X * 
a. — 
• o 

u 


— 00 . 
o — ; 

k- X 

oc ® 

-0 UJ 
X - O 
Z — ~ 

► o » 
«“ m x 
m — ea 
— a uj 
-1^-0 
Z UO X 
X > ►• 

•* < OC 
— • to cs 

o • oc 
oc 

ICOOJ 

0 ( 0 *- 


► < 
• o 


00 — 

X CN 

z - 

— cm 

OWN 

z o -« 

► Z — 

00 — (N 
O OC 4 •*• 
S ZN N 
Z • *M » 
►OX «t 
X(DX« 
ux »n 

— * — a 

X -4 Cg *- 

k- z •* o 

• XM • 

(7 *•*-(«- 

* OO -- «* 

QC tO — •» 

cs < a or 
oc a- a. — 
oc Z uj O 


*4 — UJ 

rg oo 
— -JO 
(NJ <NJ k- 
cg eg •■ 
-J X — 
X x -# 

to — * — 

•• * — to 

~ i_> O 


t-» o x 

• QC k- *- 

► -o z 

► a > - 

■ <s < o 

• K OO 00 UJ * k- k- ^ ►IN- 

► oo - z z -i •■ is — o eg » 

• uj — •* ► X — oc oo — eg ~ 

flCOOOaf (M - o »- x — 

• ► rg at uj h • -* »« — x *- 

a-zi/izrotg-iL •»—*-' 


V) Bh ►• Z * — -fNJ 


X -* • — 


-J < 

< to 


> « 
< X 
00 X 


z > — 
*“ Uj < z 

*Z W X 


00 — * 
— > < 
o qc : 

UJ S £ 

OC ►- • 


i **- Q. oO t— 


< t- CC t£ 
H 3 UJ UJ 
-SOO 
Z QC UJ UJ 
>- <C I- »- 

o z z 

oo oo — 


h- < k- 

oo x a 

► XX 

m-zz 

— I * • 

00 Q. — O 

=3 * O UJ 

aomoc 

< tm — Z 
3 k» J » 
► ►x« 
** UJ z o 
o UJ • ret 
QC w 
IL OO. 

uj m a 

~ k» 

► h w 

Z W K 
X Z -5 


<a x 
^ x 

SJ SJ 

Q O 


r z z z 
a s so si 
: x x x 

C X X X 
a o <a s> 
j u o u 


- h UJ B •» 

S&QUJ« 
a Z * oo oo 

a •- — o o 

— _j OO • M 
► ••■“> 

Z k OJ ® •• 

— O — ~ 

0 ffl - IM 

UJ Z ■ 

03 - — eg oc 

— Nl <\J ^ -» 

► ZNHUJ 

z * rg -- oO 
« Jxxa 

> X X X k- 
-J to 

UJ k Z Z Z 
X Z <3) GJ K) 

Z Z to oo to 

® ® z z z 

X X UJ UJ UJ 

1 i I I I 

B8nh- 

U U Q O O 


■ — * ►N 

-J'O'-X 
- h OJ X 
k- Si OC • 

cb oc a — 

ac • • 

-»-4 k M 

<\j — — ro 
• k- UJ -J 

•$■ es to x 

-QCOW 


x UJ UJ UJ 

o o o ►* 

ZZZZ ♦ 

a uj uj uj <■ 

- j j j g i co 

00 < < < z 

Z > > > II II It 

UJ I— — — 

•- o o a x x x 

□ 4 Ul UJ Z Z Z 


— oc — 


X -« I 

z rv eg 
l 


•t <\j r»i »r 


CO Jt 2 
X UJ UJ 

z or x. 


O 

) o z a 
: z *n z 
• — * — 
Jt UJ 2 
UJ QC UJ 


O 

z < 

— X 

□t oc 

UJ 51 
QC L. 


I uj . 

i qc s a 

- z -* z z 

I * QC — — 

* z * * 
it — <r «r 
, oc k- 

' Z to II » 
Z 

O eg eg 
O !l eg eg 


-T — -*■ 


— J JZ<f 
"5 X X k- II 
S 8 l/l l/l “} CO 
O Z — “9 Z X 


— • UJ 

OC O • 
Z UJ - 
W • OC 
k- O Z 
OO UJ — 

Z to -J 
tl Z X 
o — z 


eg — . II 
eg co — 
•J x • 
X z — 


> 

» -> < 
— » to 


* o *- -- — 

i rg — rg — 

i — eg eg oj u 

JMN k 
X -J X “> 
to ^ X > eg 
-M-< • 

»VW 1/lH 


“> to CO 
tn * a 
|| — — i ~ oc 

►- — — *4 Q. 


It ~ II 

o — — 

— o o * 

in w -4 -* 

eg —i 

O O 3 cr i 
< < -3 

a uj uj uj s * 
o ql ac x c ► 


• m — * 

i ♦ — 

> — ii 


04 ♦ — — *11 


o •— — <** uj 


' — eg co 

w Q _J KJ 

® a O eg oc 

® o I x eg co a. 
X • Z U 1 — X xz 

- O • k- OX 

— n -- n • — -- —4 

II — II — * UJ — It II 

htjCOUiZ** — kO 
» » 3 • o m 

~4 — < -4 ■— z X — eg eg 
o — o — — 

~4Q.^a.k- — ok*inm 
X X z < 

aujaujau.ujaaa 
QkOkO-XOQQ 


in in 


eo 


109 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

<u « « « o o 0 ' 0 -U'ur-r-r-r**r'-r— r-r-r'-r-<x>cocoajcococDcoaDcccrcrcrC 7 ' 

OOOOOOOOOOOOOOOCJOOOOOOOVJOOOWOOOUOO 

ooooooooooooooooooooooo«oooooooooo 



oooooooo 
^msflsooo'o^ 
q> O' O' O' O' O' o o 
oooooo-*-* 
oooooooo 


OOOOOOOOOOOOOO 0,0 o o o 

Nffl^iTO^cOO'OH CMrO ■HAOMDO'O 
O O O U O 0.0 O H H,-< M ^ r* h H (V 


oooaooooooooooo o o o o 
r- r* f*- r- r- f*- r~* r-r~ r- n- r- r*» r*~ s- r- r» r- r*- 


'»• 

1 

2 


cm in 
x m 
X '» 


CD x. — * CD 

® ui o i/> es» 

ac ►- • qc 

CL co — o » a. 

Z X I UJ -* 2 ' 


—» — ® 


-J ‘ 
* • 


“3 — • •— « H — * W 

If II — cc 

*-3 * — 

im 2 (Mcm«— . jrnin 

nk-h^z-j-t 
CM 2 1 X ~ 

^ SSSUJliSS 
X 00 Cif --«30 


—« O (| « 3 l — — 4 


j > > e» * 

» < «x o — 

:</><✓> -) 
• *• * 

> II M O — 

t « — — O i/> 

• — ■ •-* UJ z 

) <n * • < 


a/ 


m z 

O •“ 

CD 53 —I 




CO — 


II II *11 • ii — oc 2 


— X — z & * UJ 


• in 

— * 
— <r 


» • co • 


•— o z o z 


~ CM ( 
'4' 

O' UJ £ 


- (i 
OUJ -f z 
UJ • * < 

m ^ > 

« z oc 


— a w 


UJ 


C9IO si 

a <o oc 

a. m a. 

2 2 


i-ni-hZ 


•in a — «- < 

it — cc II — -« 

CD II — ' — h 2 oj |i sD II t| CT 00 2 CT 

oc>*-»-«>-cMf*“<Miii(7'— ipgajO' »z — -'•r- — — — — c^ i x — ~ •-^cn 

uj uj — «j — O' ~ o « 3 m a - cm co cm cm n j in n 

• • r— ^ r r< J z o z cm — •f^^escoiico— - j£ ® •- s 

j£ ^ « O' z iu ^ is a u; *■ s i^r-ujsNUJNUJKi; iflwznuji- 

>— — (<•. '•l-QKI-HrU-l-KZZnH^mHIC'^ tt •— - f— 

— >— * h- z — • — • •- <s *- <s» •— u. <sj — si 

u.ii'-sa.a'iaa'isssKSiau.a.siQ'ssQCScctSw'Da-'Qofo 


. -J3 — O— ^OiOCQSoO — — 0'-*003l0J •— 


CM >-* 


<M 


in in 
'T f*' 
-r <r 


in o co 

CT O* 


•>4 in N <0 cci o 

r- IwNNI*. CO 

rn i« cn i*i n cn 


II do — 

“^oinzo 1 

? CM 3 

n •— Z Q^ 

OD UJ « Z O 
mt-K-NQ 
' < 

S) 1 S UI S UI 
Q1UCC3CC 

U1KI- O' 

CD CO O' O' 

ro cn m o' 


— . j • ,i 

M 2 £* 

O II. 2 C"** 


110 



\ 


• OOOOOoOO O O O © ©'© 00000.00 o o oooooooooooooooooooooooooooooooo o'o o o o o 
^r\jf0^u>^)r i ^oDO'0^cM<^>fin-cr-aDcro»-<rsiro<'irt'Cr»cocro<-<f\jpo-4-ir>>cr~cDcro-*rNjn*4-m>or-cDcro-Hfvjc^>riri'Or-aoo'0*M 
, njnjnjfgc>jcgr^nj^rnrnrnrn(^oirnfrifnrn^»r<^^<r'T^'r^tni/>ininininin»nin'inv©»o>0'0^gj«o«o^J'Or“P' , *r>>-r»r“r>“r-n“f'*r , *-aocc 

O O'O O O O O OO O O O O O O O © O OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 





X 

« 




® “5 wizzo * a * * x 

O * V) o v O !✓» • CO i— CD < CD —03 <— ' 4 ' </>!✓) OO »CD 

*-• • o i i i hs -a uj® -a - •■ u u z z • s> -• ea 

— . — - 4 - >r -r ii h it it ii a " 4 - ^t/xt'roaa^ co uj oc lu -*■ *■ **■ I i^whoc 4 - 4 - *— at «* 


*- "• * o o r» r- -^a)uj xxoxo--**- — »«-»a.xoxoax* aioxviaiMizozoiiniiii Haxoihax 

i/> ^ <n -• *- in in • ..go zz • z •ra^^rgmzz • z*~zzo-*zz • z o z z o z z • z z z • z s z 2 

ZZUMhl- •• O w o «»0 •>•-‘44444 »»0»Q»» O *> ► O ► ► H •• •> O » -* <\J ^ r\J H"! * - O •■ 0C •• *■ 

♦ ~4 •-* Ol — -* f-t II II' “5 “5 “5 “) -J ^4 ^ It — i II II I 44444^-* 

uo»-hz w • m k 11 — n —. •■ ► *> » •> ii 11 11 — 11 n , o rn ~ 11 - 11 11 ~ n 11 — u n "s “> -5 -> -7 11 11 h —• 11 11 

UMUXVNai/M/lZIT'OOUjO - * * » • 

i/im«>—i/}oc 4 ®>— «\i*»ouj*-ogo • *44444 » » —) • * « » * th ^ „ (\j m * » 

«*-•»/> 00 O l/> ^ o z Ifl -ZNtMwIfmmTJO'C-'O — O O »X(MN~Nw|f -4 4 4 4 4 U1 in •— lA — t\f f\J 

H H 5 J IS ® -I Z CO «- 00 >— O-w.-,— — O — irvin-ZOO— O'-OO'-OO— ' 0 " 5 ">“ 5->“>->00 — o~oo 
O O It »— II || H It II H II IIK»NO|-m(nOinOQaOOQWf , tW^^f , t' , l*'-f'IMlUMU;NNliJinWQ(M*'>-'- w *.wN(Ml/lNl/)'f'f 
W UJ O z < Z h- ►- »- h* H- h- OO *- I/O lO in 4 - H- K- ^ H- K ►— OO 

</)i/iNisizi/)SZi/)(az^^u.siujiais(ssQsiasiaiss) 5 }ca M 5 ) M SS)jLLcasQsos 6 ios> 5 ississtsisssosM 5 )Msis 
"“l-lS^OOiflUOWUO-QaUOOKQatfl'aaiKXDaDCOQOX-QCt-OI-CQl-ClCXQacfXj'arSOD'i.CU.CO 


CO #4 ■’■' IN ft ■ Ifl 

cm • ry fM n# 


co © r\i • * ^ in 

r- co o o 

m fr\ 


■C O ' nj p r 

©in o o 

cn m cm rg 


«m in - 

O O > O 

m pg f\i 


111 






11 ? 



SUBROUTINE LEBEGE 


Hie subroutine link LEBEGE receives the YICS array for each segment from 
[NITAL via Tape #4. The subroutine FIXEM is called to integrate the 
differential equations of each segment, under true load conditions. 

FIXEM is identical to subroutine SETUP, while WAND corresponds to subroutine 
MAGIC and only consideration of the OVERLAY structure dictates the change in 
names. The subroutines TOBAR, TEMQEG, PLYCO, and PLYNE are similarly 
equivalent to ROBOT, GEOMET, PLICO, and PLINE discussed previously. 

The results of the final integration sequence are the forces and 
deflections at the beginning, intermediate, and end points of each segment. 
In the first pass, first cycle, the prestress and predeformation states are 
printed from this routine. In other passes the U, V, W eigenvectors are 
printed, providing convergence has been reached. 
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SUBROUTINES ODE1 AND 0DE2 


Subroutine LEBBGE calls either ODEl or ODE2, as necessary, and various geometric 
and trigonometric clues, as well as the predicted values of the variables for the 
differential equations, are passed to this subprogram via label common area LASTEQ. 

The equations in ODEl and 0DE2 are identical to those in subroutines DIF1 and 
DIEF2 respectively, with the addition of the auxiliary equations for YAQPH, and 
YAQTH. Subroutines ODEl and 0DE2 perform the final integration for each segment 
in the structure utilizing the initial conditions previously obtained, and return 
these values to LEBEGE via label common area LASTEQ. 





The ODEl, 0DE2 flow charts are identical to the DIF1, DIFF2 flow charts, respectively 
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SUBROUTINE EIGVAL 


This is the controlling routine for the program eigenvalue calculation 
loop. The necessary matrices for the calculations are passed to COMPAK, 
and thence to routine EIGEN. The controlling routine EIGVAL will 
provide intermediate matrix print as desired. 

Subroutine COMPAK : This routine reads rows of a two-dimensional array 

from a storage unit into core and then stores them into a column vector. 
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SUBROUTINE EIGEN 


EIGEN is the main eigenvalue and eigenvector calculation routine. 

Printout of the eigenvalues and appropriate eigenvectors is provided by 
this routine. First EIGEN prepares the matrices for eigenvalue extraction, 
calling FUTILE, DAGGER, and SWITCH. Then the matrices are passed to 
SYMEIG, which obtains the eigenvalues and specified vectors using the 
Householder technique. 

Subroutine DAGGER ; This routine performs the reduction of the two matrix 
eigenvalue problem to a standard eigenvalue problem. 

Subroutine SYMEIG (Alternate entry point SYMVEC) : This routine obtains 

the eigenvalues and vectors of a real, symmetric, matrix. . 

Subroutine TFORM: This routine performs the Householder similarity 

reduction of a real, symmetric matrix to tri-diagonal form. 

Subroutine STURM: Performs the Sturm-sequence calculations of eigenvalues 

of a real, symmetric, tri-diagonal matrix. 

Subroutine PREP (alternate entry point PET) : Determines the number of roots 
greater or equal to a stipulated value. 

Subroutine QSVFC (alternate entry point QWIEL) : Obtains eigenvectors of a 

tri-diagonal matrix and back transforms them to eigenvectors of the original 

1 

real, symmetric matrix. 

Subroutine RANDOM : Random number generator. This subroutine is an MSFC 

system subroutine and the listing is not included here. Most any random 
number generator may be substituted or the argument may be set to zero. 

Subroutine ANDD: This routine "ands" or "ors" the corresponding bits of 

two given variables with each other. 

Subroutine DOTPRO : Performs inner product accumulation using labeled 

common area INFO. 
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50 CONTINUE _ 3600840 

REWIND 4 • { ' 3600841 

99 RETURN . ^ ,3600850 

end 3600860 
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SUBROUTINE DETERM 


This is the controlling routine for the program determinant evaluation 
loop. The matrix in question is passed to DET2 where the determinant 
evaluation calculations occur. 

Subroutine DET2 : This routine computes the determinant of a square matrix. 

The actual determinant value is normalized to ±1. 

This routine performs the triangular factorization of 


This routine performs inner product accumulation. 


Subroutine DC0MP2: 
a square matrix. 

Subroutine SUPER: 



CALL DCOMP 2: PERFORM TRIANGULAR 
FACTORIZATION OF A SQUARE MATRIX- 
CALLS SUPER TO ACCUMULATE THE 
INNER PRODUCT 
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SUBROUTINE BCVECT 


This routine converts the eigenvector obtained during a converged 
eigenvalue pass into the form accepted by INITAL. The calculations 
are identical to those in the latter part of STRMAT, and result in 
the DRE array. 
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BCVECT 
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SUBROUTINE ETRAP 


This is an error trap routine which can be called by the MAIN routine 
at various stages of program execution. If the indicator NIX is not 
equal to zero, MAIN will call ETRAP and indicate the proper message 
to be printed. 
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